The Magnetic Structure of DyMn 2 05 Determined by Resonant X-ray Scattering 
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Resonant magnetic x-ray scattering has been used to investigate the magnetic structure of the 
magnetoelectric multiferroic DyM^Os. We have studied the magnetic structure in the ferroelectric 
phase of this material, which displays the strongest ferroelectric polarisation and magnetodielectric 
effect of the RM^Os (where R is a rare earth ion, Y or Bi) family. The magnetic structure 
observed is similar to that of the other members of the series, but differs in the direction of the 
ordered moments. In DyM^Os both the Dy and Mn moments lie close to the b axis, whereas in 
other RM^Os they lie close to the a axis. 
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I. INTRODUCTION 

Magnetoelectric multiferroics have been widely stud- 
ied for the last half century because of the interesting 
physics they display^—, and the prospects for applica- 
tions based on the mutual control of magnetic order and 
ferroelectricitjsi. Much of the recent focus has been on 
the properties of TbMnOy^ and similar materials with Tb 
replaced by another rare earth ion 7 -, and on the RM^Os 
series (where R is a rare earth, Y or Bi). The mecha- 
nism behind the multifcrroicity in these systems is not 
completely understood, but theories based on symmetry 
considerations have provided important insights^. To 
construct microscopic models that apply to these sys- 
tems more information is needed on the details of the 
various magnetic and electronic phases. 

The magnetic properties of many members of the 
RMn2 05 series are well documented through various 
studies going as far back as the 1960 o 10 ' n . There has 
been a recent resurgence of interest in RMn205 com- 
pounds since the discovery of the magnetically-induced 
ferroelectricity in these materials^ 2 - - — . The specific 
magnetic structures are believed to be important to 
the magnetoelectric coupling and have been established 
for many members of the series, usually by neutron 
diffractio n 1 1 ' 17 i l2. or by resonant x-ray scattering 2 ^. 

Among the RMi^Os series, DyM^Os displays 
the largest electric polarization 2 ^ and magnetodielec- 
tric effect^. Previous studies of the magnetic 
structur e 11 ! 17 ' 22 ^ 3 have examined the low temperature 
phase and followed the magnetic ordering wavevector as 
a function of temperature, but crucially, the magnetic 
structure in the ferroelectric phase has not been fully de- 
termined. Here we address this outstanding issue by re- 
porting the magnetic structure of DyM^Os in the ferro- 
electric phase as determined by polarized resonant mag- 
netic x-ray scattering (RMXS) and ab initio calculations 
of the resonance spectrum. 

The bulk properties of many of the RM^Os com- 
pounds are well documented. They show a ferroelectric 



polarization along the &-axis, approximately in the tem- 
perature range 20 < T < 35 K. There exists a transi- 
tion Tjv ~ 40 K into an incommensurate magnetic phase 
(HT-ICM), just above the ferroelectric phase. Below the 
HT-ICM phase, the magnetic structure becomes com- 
mensurate and is coexistent with the onset of ferroelec- 
tricity (CM-FE phase). A low temperature incommensu- 
rate phase (LT-ICM) appears below 20 K. These phases 
only have spontaneous magnetic order on the Mn ions, 
but at lower temperatures still (T < 10 K) another phase 
exists having spontaneous ordering of the magnetic rare 
earth ions. 

In DyMn 2 5 the HT-ICM phase occurs between 42 
and 44 K, and the CM-FE phase exists down to ~15K. 
The LT-ICM phase is observed below 20 K. DyMn 2 5 
displays some differences from the other members of 
the RMii2 05 series. In addition to small variations in 
the phase transition temperatures, the ordering wavevec- 
tors of the incommensurate phases are different, and 
DyM^Os is unique in having a commensurate structure 
for the spontaneous order of the rare-earth moments* 2 ^. 

DyMn2 05 has a room temperature orthorhombic 
structure, space group Pbam, with lattice parameters 
a = 7.294 A, b = 8.551 A and c = 5.688 Ail. The Mn and 
O form a structure of Mn 4+ 06 octahedra and Mn 3+ 05 
square-based pyramids, with the Dy ions located in the 
gaps of the MnO structure, as shown in Fig.[TJ The mag- 
netic structure of DyM^Os is fully established below 
8 K 11 ' 17 . Neutron diffraction 2 ^ and RMXS 2 - experiments 
were made in the CM-FE phase and confirmed that the 
magnetic scattering occurs at (h + 1/2, k, I + 1/4) reflec- 
tions, but these studies did not determine the magnetic 
structure in this phase. 

The large neutron absorption cross-section of natural 
Dy limits what information can be obtained by neutron 
diffraction on the magnetic structure of DyM^OV As 
an alternative to neutron diffraction, RMXS has two im- 
portant strengths. Firstly, the resonant enhancement in 
the scattering when the photon energy is tuned to an 
atomic transition enhances the otherwise very weak x- 
ray scattering signal from the magnetic order, and sec- 
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FIG. 1. The crystal structure and magnetic structure (model 
A) of DyM^Os viewed along two different crystallographic 
axes, (a) shows the structure when looking at it along the 
c-axis and (b) looking at the structure along the b-axis. In 
both images, the Dy ions are the blue circles, the Mn 4+ ions 
are magenta and are positioned inside the oxygen octahe- 
dra and the Mn 3+ are shown in green and are positioned 
inside the oxygen square-based pyramids. The unit cell is 
shown by the black box. In (b) only the moments on the 
Mn 4+ sites are shown, demonstrating the interaction between 
two Mn 4+ moments along the c-axis within a unit cell is fer- 
romagnetic. 



ond, a degree of sensitivity to the separate order on the 
Dy and Mn sublattices can be obtained from the vari- 
ation in scattering at different absorption edges. These 
strengths were exploited in the recent RMXS study of 
magnetism in DyMn03. 24 



II. EXPERIMENTAL PROCEDURE 

RMXS was performed at the ESRF, Grenoble. The 
beamline ID2QS5., a high-intensity beamline with full po- 
larization analysis, was used with a four-circle diffrac- 
tometer in vertical scattering geometry, shown in Fig. [3] 
A high quality single crystal of DyM^Os, grown in Ox- 
ford et al^, was used for this experiment. The crystal 
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FIG. 2. The vertical scattering geometry used in this experi- 
ment. 



was cut so that the (0,0,1) direction was perpendicu- 
lar to the largest face of the crystal. The crystal was 
mounted on the cold finger of a closed-cycle refrigerator 
to allow cooling as low as 12 K. 

Measurements were performed at two x-ray energies, 
corresponding to two resonance energies in the ions be- 
ing studied. These were the Dy Lni-edge (2p — > 5d), and 
the Mn isT-edge (Is — ¥ Ap). The measurements were per- 
formed at the peak intensity of the resonance spectrum. 
Fig. [3] shows the energy line shape in the vicinity of the 
Dy .Lin-edge, measured in the an' polarization channel 
(which shows a purely magnetic signal) at a tempera- 
ture of 15 K, with the x-ray scattering vector fixed at the 
(—0.5, 0, 4.25) magnetic Bragg peak. The inset to Fig. [3] 
shows the fluorescence signal in the same energy range. 

Polarization analysis was used to determine the mag- 
netic structure, following the theoretical descriptions of 
Blume et al^L and Hill et al£&. The intensity and polar- 
ization of the resonant x-ray diffraction signal is depen- 
dent on the magnetic structure. Two different methods of 
analysing the polarization dependence of the x-ray scat- 
tering from the magnetic structure were employed. These 
are discussed in the following paragraphs. 

Azimuthal scans of certain magnetic Bragg reflections 
were measured. In this type of scan the incoming lin- 
ear polarization is fixed perpendicular to the scattering 
plane (tr polarisation) and the crystal is rotated around 
the scattering vector. The analyzer is set so that only 
scattered x-rays with polarization within the scattering 
plane (tt' polarization) are detected. The integrated in- 
tensity (measured using a 0-scan) of the magnetic Bragg 
peak in this <77r' channel is then obtained as a function 
of azimuthal angle. This technique has been used to de- 
termine details of the magnetic structure in the related 
materials TbMnO^^, HoMnzOg^ and TbMn 2 0g2&. 

All of the measurements on the CM-FE were con- 
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FIG. 3. A scan of the energy dependence of the (-0.5,0,4.25) 
magnetic peak, through the resonant Dy Lm-edge, at 15K in 
the an' polarisation channel. The measured data is shown 
with black circles and the line is the energy profile for model 
A calculated by FDMNES. The inset shows the fluorescence 
measurement around the same edge, again with the measured 
data shown with black circles and the FDMNES calculation 
shown with a black line. 



ducted at 15 K. In total 5 azimuthal scans were measured, 
four at the Dy Lm-edge, at wavevectors (—0.5,0,4.25), 
(-1.5, 0, 3.25), (0.5, 1, 4.25) and (0.5, 2, 5.25), and one at 
the Mn K-edge, at wavevector (0.5,1,3.25). For all re- 
flections, the azimuth reference vector defining the zero 
of the azimuth was [0,0,1]. The Dy Pm-edge scans are 
shown in Fig. |4]and[5j an d the Mn LT-edge scan is shown 
in Fig. [6] together with an energy scan through the weak 
K-edge resonance (insert to Fig. EJ)- The data have been 
corrected for the absorption of the x-ray beam in the sam- 
ple, which varies with azimuthal angle. 31 This correction 
is essential because the absorption varies strongly with 
the angle that the incoming and outgoing x-ray beams 
make with the crystal surface. 

The alternative method is to use full linear polariza- 
tion analysis (FLPA), in which the incoming x-ray po- 
larization is varied by use of phase-plates, while keep- 
ing the sample in a fixed position ! 20 ' 32 This has two ad- 
vantages over the azimuthal scan. First the attenuation 
does not vary during the measurement, and second, the 
sample does not move so the same part of the sample is 
measured. The sample was aligned on the (—0.5,0,4.25) 
reflection at 15 K, again using the resonance at the Dy 
Lni-edge. The azimuthal angle the scan was performed 
at was —97.5°, the maximum of the intensity in the az- 
imuthal scan at this wavevector. The angle of the in- 
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FIG. 4. Azimuthal scans performed at the Dy Lm-edge, at 
15K, in the on' polarisation channel. The top panel shows a 
scan at the (—0.5, 0, 4.25) magnetic reflection, and the bottom 
panel shows the ( — 1.5,0,3.25) magnetic reflection. The data 
points are shown with filled circles and the fit to magnetic 
structure model A is shown by the line. 



coming polarization, ?y, is set so that r\ = corresponds 
to a polarised light (and therefore tt polarised light is 
at rj = 90°). For each value of incoming polarization 
measured, a FLPA of the outgoing beam was performed 
to obtain the values of the first two Stokes parameters, 
Pi and Pi. The parameter Pi describes the polarization 
perpendicular and parallel to the scattering plane and P2 
the polarization at ±45° to the scattering plane. The two 
parameters th erefore ch aracterise the linear polarization, 
with P Un = y/P? + P*. There is a third parameter, P 3 , 
which is related to the circular polarization. Fig. [7] shows 
Pi and Pi obtained this way as a function of 77. 



III. RESULTS 

To determine the magnetic structure of DyM^Os, we 
treated the analysis in three parts. The magnetic order- 
ing on the Dy 3+ and Mn 4+ sites was determined from a 
fit to the four azimuthal scans and the FLPA 77 scan per- 
formed at the Dy Lni edge. Having obtained the Mn 4+ 
magnetic structure, we determined the ordering on the 
Mn 3+ ions from a fit to the azimuthal scan at the Mn K 
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edge. Finally, we confirmed the validity of the results by 
comparing the energy scans over the Dy Lni edge with 
ab initio calculations. The procedure is explained in the 
following paragraphs. 

X-ray magnetic scattering at the Dy Lm edge probes 
the empty Dy 5d states. As has been discussed 
previously^ the magnetic polarization of the 5c? states 
arises not only from the localized Dy 4/ moments but 
also from the Mn 4+ moments due to the overlap between 
the extended Dy 5d states and the Mn 4+ 3d orbitals. The 
Mn 3+ ions are thought to have a negligible contribution 
at the Dy Lm edge, as they are situated further from the 
Dy ions in the crystal structure. 

Previous work has shown that at low temperatures the 
spontaneous magnetic order of the rare earth sublattice 
does not always have the same symmetry as that of the 
Mn sublattices2&22. Following work on TbMn 2 5 and 
HoM^O j 18 ' 20 , we constrained the direction of the Dy 3+ 
moments on pairs of sites to be the same. Specifically, 
sites (1,2) have one direction and sites (3,4) have another 
direction (see Table Q] for the positions of the sites), and 
we constrained the size of all the Dy 3+ moments to be 
the same. For the Mn magnetic order we assumed the 
same structure as found many times previously 49 . This 
assumption has the Mn 4+ moments linked ferromagnet- 
ically along the c-axis within a unit cell. The Mn 4+ ions 
at x — and x — 0.5 are linked by a &-glide. The Mn 3+ 
structure has two Mn ions linked by a 6-glide and the 
other two linked by a fo'-glide. Initially we neglected the 
component of the ordered moments along the c direction, 
which is reasonable given that the magnetic structures of 
other RMn205 compounds have only a small c-axis com- 
ponent. 

The fit to the data at the Dy Lm-edge was performed 
simultaneously to the four azimuthal scans and the FLPA 
j] scan (Figs. [H [5] and [7|). The free parameters were 
the three in-plane angles of the ordered moments (one 
for each of the inequivalent Dy sites, and one for the 
Mn 4+ site) and the fraction a of the Mn 4+ contribution to 
the signal (the remainder coming from the Dy 3+ ). The 
best fit was achieved with a = 0.64(7). We note that it is 
not possible to obtain a good fit to the Dy Lm-edge data 
allowing only Dy 3+ ions in the model (a = 0), nor is it 
possible using only Mn 4+ ions in the model (a = 1). This 
demonstrates that the signal at the Dy Lni-edge arises 
from both Dy and Mn magnetism. 

At the Mn K-edge, the empty Ap states of the Mn 
ions are probed, but the magnetism is carried by the 3d 
electrons. The Mn K-edge scan is sensitive to the Mn 4+ 
and Mn 3+ moments via an intra-atomic Coulomb inter- 
action. We assume the contribution from the Dy mo- 
ments is negligible because the Dy 4/ and Mn 4p states 
are highly localized. We fitted the Mn K-edge azimuthal 
scan to a model for the magnetic ordering of the Mn 3+ 
ions assuming the magnetic structure of the Mn 4+ ions 
determined previously from the data at the Dy Lni-edge. 
In principle, it should be possible to distinguish between 
Mn 3+ and Mn 4+ signals using the characteristic energy 
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FIG. 5. Azimuthal scans performed at the Dy Lm-edge, at 
15K, in the an' polarization channel. The top panel shows a 
scan at the (0.5, 1,4.25) magnetic reflection, and the bottom 
panel shows the (0.5, 2, 5.25) magnetic reflection. The data 
points are shown with filled circles and the fit to magnetic 
structure model A is shown by the line. 



line shape of the Mn K-edge resonance. Unfortunately, 
the K-edge enhancement is very weak, and the data are 
not good enough for this type of analysis (Fig. [6] insert). 
Instead, we determined the ratio /3 of the contribution 
from Mn 3+ to Mn 4+ by the same method as described 
in the previous paragraph for a. The best fit is shown in 
Fig. [6]and had /3 = 0.55(2). 

Using this method, several equally good fits of the az- 
imuthal and FLPA T) scans were found. However, all 
but one of the models could be dismissed from analy- 
sis of Dy Lin-edge resonant energy spectrum. To dif- 
ferentiate between these magnetic structure models, ab- 
initio calculations of the Dy Lni-edge resonance at the 
(-0.5,0,4.25) reflection were performed. The FDMNES 
code was employed 3 ^, which performs a monoelectronic, 
cluster-based multiple scattering calculation, well suited 
to the prediction of resonant spectra that involve excita- 
tions into delocalized (rare-earth 5d) states. A 256 atom, 
2a x b x 4c supercell (the magnetic unit cell) was de- 
fined, containing 32 Dy atoms, located on 8 symmetry 
inequivalent sites. A cluster radius of 4.5 A was selected, 
above which no improvement in the calculation accuracy 
was observed. First, the absolute energy scale was cali- 
brated by calculating the Dy Lm absorption spectra and 
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FIG. 6. Azimuthal scan performed at the Mn Tf-edge, at 
15K, in the an' polarization channel. The wavevector is 
(0.5,1,3.25). The experimental data are shown with filled 
circles, and the fit to magnetic structure model A is shown 
by the line. The insert shows the energy line shape of the Mn 
K-edge resonance at the same wavevector and temperature, 
measured at an azimuthal angle of -79°. The resonance peak 
is marked by the black arrow. 



comparing to the experimentally measured fluorescence, 
shown in the inset of Fig. [3] The Fermi energy was then 
set at Ep = 7.778 keV, giving a low energy cut-off. The 
resonance was calculated with both El-El and E2-E2 
excitation channels allowed. The E2-E2 excitation, how- 
ever, gave a negligible contribution to the resonant in- 
tensity and was therefore omitted from the final calcu- 
lations. By comparing calculations based on the possi- 
ble magnetic structure models it was clear that only one 
model comes close to reproducing the resonant energy 
spectrum, which is shown in Fig. [3] 

The magnetic structure which is consistent with the 
azimuth scans, FLPA scans, and resonance lineshape is 
presented in tableland depicted in Fig.[TJ This model of 
the in-plane magnetic structure will hereafter be referred 
to as model A. As mentioned earlier, we only fitted the 
in-plane components of the magnetic moments. The di- 
rections of the moments are therefore specified by their 
angle 4> from the a-axis. The layers stack with alternat- 
ing ferromagnetic and antiferromagnetic coupling along 
the c axis. This gives the ttJ-4- propagation of the mag- 
netic structure along the c-axis that is observed in other 
RMn2C>5 materials. 

The fact that the temperature at which these exper- 
iments were performed (15 K) is above the spontaneous 
ordering temperature for Dy (8 K) suggests that any or- 
dered magnetic moments on the Dy sites are induced by 
the magnetic order of the Mn sublattice. In this scenario, 
it is expected that the ordering on the Dy sites should 
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FIG. 7. A full linear polarization analysis measurement of 
the Stokes parameters Pi and P2 as a function of incoming 
polarization (rj). The data were collected at the Dy Lm-edge 
at a temperature of 15 K. The X-ray scattering wavevector is 
fixed at (—0.5, 0, 4.25) with the angle of the azimuth at -97.5°, 
the maximum of the azimuthal scan at this wavevector. The 
symbols represent the experimental data, and the lines are 
calculated from the magnetic structure model A. 



follow the symmetry of the Mn magnetic structure. How- 
ever, this symmetry constraint is not applied in model A. 
Therefore, we tested a second model (model B) with the 
moments on the Dy ions following the same symmetry as 
the Mn moments. 

In model B, the Dy moments are still considered in 
pairs, but rather than being parallel, as in model A, the 
moments in each pair are related by either a 6-glide or 
fe'-glide. Under these constraints the best fit to the data 
at the Dy Lni edge is found when a = 0, meaning that 
only the magnetic moments at the Dy sites contribute to 
the signal. The calculated azimuthal and FLPA 77 scans 
for this fit are very similar to the ones presented in fig- 
ures HI [5] and [3 and are therefore not shown here. The 
value of the reduced x 2 is about 3% higher for the best 
fit to model B when compared to the best fit to model 
A. An FDMNES simulation (not shown here) assuming 
the best-fit model B structure is consistent with the mea- 
sured Dy Liu edge resonance at the (—0.5, 0, 4.25) reflec- 
tion. Fig. [5] depicts the model B magnetic structure. 

The experimental geometry used meant that we had 
limited sensitivity to the component of the magnetic mo- 
ment parallel to the c-axis. Test models in which a c-axis 
component was added by hand confirmed this insensitiv- 
ity, which is caused by the relatively small angles that 
the scattering vectors at each of the reflections studied 
make with the c-axis. We cannot, however, rule out a 



6 



TABLE I. The magnetic structure (model A) of DyMn 2 5 at 
15 K, in the commensurate antiferromagnetic and ferroelec- 
tric phase. The positions of the ions have been determined 
elsewhere^. The measured moment directions lie in the aft- 
plane and are specified by the angle (j> from the a-axis. The 
numbers in brackets in the final column are the uncertainties 
in the last digits obtained from the fits. 
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FIG. 8. Magnetic structure model B of DyMn205, in which 
the Dy magnetic structure is constrained to have the same 
symmetry as the Mn magnetic structure. 



small c-axis component to the ordered moments, but as 
noted earlier, in other RMn2 05 compounds the moments 
lie predominantly in the a6-planc. 



IV. DISCUSSION 

The magnetic structure of the Mn ions, found here 
for the ferroelectric phase, is similar to the magnetic 
structure of the Mn ions found in DyM^Os in the Dy- 
ordered phase by previous studie s 11 ' 17 . The structure 



shows approximately antiferromagnetic chains running 
parallel to the a-axis, a common feature of the RM^Os 
series that agrees with the magnetic exchanges discussed 
elsewher o 17 ' 35 . However, one distinct difference between 
DyM^Os and the other members of the RM^Os se- 
ries is that in the ferroelectric phase measured here, the 
magnetic moments are aligned close to parallel with the 
6-axis. This is in contrast to the other members of the 
series which tend to have the magnetic moments approx- 
imately parallel to the a-axis. 

The best fit model we obtained (model A) breaks the 
space group symmetry on the Dy site. In this model, 
the Dy moments occur in ferromagnetically-aligned pairs. 
This is similar to the rare-earth magnetic ordering in 
TbMn 2 5 ^S, but in contrast to that in HoMn 2 3 1 ^. 
However, at 15 K, the material is above the spontaneous 
magnetic ordering transition temperature of the Dy ions, 
so at this temperature it would be expected that any 
static order of the Dy moments would be induced by 
coupling to the Mn sublatticc and would therefore have 
the same symmetry as the Mn magnetic order. This is 
the case for model B (by construction) but not for model 
A. Therefore, if model A describes the correct structure, 
then the implication is that the Dy-Dy magnetic inter- 
actions remain a sufficiently strong influence at 15 K to 
prevent the Dy magnetic order from adopting the same 
symmetry as the Mn magnetic order. 

The differences between the Dy magnetic structures of 
model A and model B are actually relatively small. In 
model A, the Dy moment on site 4 is parallel to that 
on site 3, whereas in model B the moments on sites 3 
and 4 are approximately antiparallel (see Figs. da) 
and [8]). The other Dy moments are found to have only 
a slight difference in their angle to the a-axis between 
the two models. In both models the Dy moments are 
close to parallel with the 6-axis, again in contrast to the 
observations of other members of the RM^Os series. 

There are a few noteworthy features of the fit itself. 
The fit successfully describes three out of four of the Dy 
Lni-edge azimuthal scans, but does not adequately de- 
scribe the azimuthal scan at (0.5, 2, 5.25). When treating 
this scan independently of all other scans, we could find 
no model that describes the observed azimuthal depen- 
dence. One possibility is that there is a systematic error 
due to sample absorption. However, calculations indi- 
cate that the variation in absorption is not significantly 
greater than at the other reflections. 

In the FLPA ?7-scan, around the point r\ = 90°, Pf 
and P| both tend to zero. The magnitude of the Stokes 
vector (Pi, P2, P3) should be equal to one for polarized 
light. It is obvious that at 77 = 90°, P Un ^ 1 (Fig. 0). 
The data suggests that the beam is either depolarized or 
there is a circular component to the polarization. This 
effect has been observed before and was explained by a 
non-zero P3 and therefore a circular component to the 
polarization of the outgoing bear m 32 ' 36 ' 37 
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V. CONCLUSION 



In this work we have investigated the magnetic struc- 
ture in the multifcrroic phase of DyMnaOs, which has 
not been considered in detail before now. We have found 
two alternative but similar models which provide a good 
description of the data. Although the symmetry of the 
rare-earth and Mn magnetic structures in DyM^Os is 
the same as that found in other RM112O5 compounds, a 
key difference is that in DyM^Os the Dy and Mn mo- 
ments point in a direction close to the b axis, whereas 
in other RMn 2 05 compounds the moments are close to 



the a axis. Given the strong coupling between magnetism 
and ferroelectricity in the RMn 2 05 series, it is likely that 
this difference plays a role in causing the exceptionally 
large electric polarisation in DyMi^Ov 
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